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Abstract 

We analyze a class of dynamics of open quantum systems which is governed by the 
dynamical map mutually commuting at different times. Such evolution may be effectively 
described via spectral analysis of the corresponding time dependent generators. We consider 
both Markovian and non-Markovian cases. 

1 Introduction 

The dynamics of open quantum systems attracts nowadays increasing attention [IJ. It is very 
much connected to the growing interest in controlling quantum systems and applications in 
modern quantum technologies such as quantum communication, cryptography and computation 
[2j. The most popular approach is to use a Markovian approximation and to consider a master 
equation 

At, to = ^tAt,to , ^to.to = id , (1.1) 

with time dependent generator Ct- The above equation gives rise to a quantum dynamical map 
(completely positive and trace preserving) At which in turn produces the evolution of a quantum 
state pt = Atp. The corresponding generator Ct has to satisfy well known condition [3', ^ (see 



also i5j for the detail presentation) and the solution is given by the following formula 



=Texp^ Cudu^ , (1.2) 

where T denotes a chronological product. We stress that the above formula has only a formal 
character since the evaluation of its r.h.s. is in general not feasible. If the generator does not 
depend on time Lt = C then it simplifies to 

At^to = exp (£ (t - to)) • (1.3) 

Let us note that characteristic feature of ()1.3p is that Markovian semigroup At := At-\-to,to is 
commutative, that is 

AtA, = AsAt, s,t>0. (1.4) 



It is no longer true for the general time dependent case (II. 2p . The general Markovian evolution 
does satisfy the inhomogeneous composition law 

-^t,u-^u,s — -^t,S 1 (1-^) 

for t > u > s > to, however, it is in general noncommutative. 

Non-Markovian evolution is much more difficult to analyze (see [6]-|17] for the recent papers). 
The local master equation is replaced by the following equation 

At,to = / K^t-uAu,to du , p{to) = po , (1.6) 

Jto 

in which quantum memory effects are taken into account through the introduction of the memory 
kernel /Cj: this simply means that the rate of change of the state p{t) at time t depends on its 
history (starting at t = Iq). Recently, we proposed a different approach [18] which replaces the 
non-local equation (II. 6p by the following local in time master equation 

Atfy = Ct^tAtM , ^to,to = id • (1-7) 

The price one pays for the local approach is that the corresponding generator keeps the memory 
about the starting point 'to'- This is the very essence of non-Mar kovianity. Interestingly, this 
generator might be highly singular, nevertheless, the corresponding dynamics is perfectly regular. 
Remarkably, singularities of generator may lead to interesting physical phenomena like revival 
of coherence or sudden death and revival of entanglement [18]. Now, the formal solution to (II. 7p 
reads as follows ^ ^ 

At^t, =Texp( / ' Cudu] . (1.8) 



It resembles very much Markovian dynamical map ()1.2p and again its r.h.s. has only formal 
character due to the presence of the chronological operator. Note, however, important difference 
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between ()1.2I) and (II. 2^ : the former does satisfy composition law. The latter is homogeneous in 
time (depends upon the difference t — to) but does not satisfy ()4.26p . 

In the present paper we analyze a special case of commutative dynamics, i.e. we generalize 
(|1.4|) for time dependent Markovian and non-Mar kovian dynamics. In this case formulae (|1.2p 
and (|1.2p considerably simplify - the chronological product drops out and may compute the 
formula for the dynamical map via spectral analysis. 



2 Preliminaries 

Consider d-dimensional complex Hilbert space and let {eo, • • • ,6^-1} be a fixed orthonormal 
basis. For any x,y ^ denote by {x, y) the corresponding scalar product of x and y. Let 
Md = C{C^, C^) denote a space of linear operators in C^. Now, is equipped with the 
Hilbert- Schmidt scalar product 

d-l 



(a, b) := ^(aefe, bet) = tr(a*6) , (2.1 
where a* : C — t- is defined by 



k=0 

is defined by 

{a*x,y) = {x,ay) , (2.2) 

for arbitrary x,y ^ C^. Finally, let us introduce the space C{Md, M^) of linear maps A : — t- 
Md- For any A £ C{Md, M^) one defines a dual map A#eC{Md,Md) by 

{A*a,b) = {a,Ab) , (2.3) 

for arbitrary a,b ^ M^. Note, that if the dual map A"^ is unital, i.e. A^I^ = Id, then A is trace 
preserving. It is clear that C{M(i, M^) defines d? x d? complex Hilbert space equipped with the 
following inner product 

d2_l d2_i 

{{A,B)) = Y,{Af^,BU)= Y^trUUnBf^)] , (2.4) 

for any A,B& C{Md, M^)- In the above formula fa denote an orthonormal basis in M^. Let us 
observe that in C{Mci, M^) one constructs two natural orthonornal basis 

: Md^Md, (2.5) 

and 

: Md^Md, (2.6) 

defined as follows 

Faf^a = faafp , (2.7) 

and 

Eapa = faUp.a) , (2.8) 
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for any a G M^. One easily proves 

{{Fa^, F^,)) = {{Ea^, E^,)) = 5ai,5fS. ■ (2.9) 
Moreover, the following relations are satisfied 

^ F^aa = Idtra , (2.10) 

a=0 

and 



^Ec,aa = a. (2.11) 



a=0 

Remark 1 Note, that representing a linear map A in the basis Fq,^ 

A = J2aapFa^, (2.12) 

with 

aap = {{A,Fa^)) , (2.13) 

one has a simple criterion for complete positivity of A: a map A is complete positive if and 
only if the corresponding cf x matrix ||aa/3|| is semipositive definite. On the other hand the 
^^-representation 

with 

a'^^ = {{A,E^fs)) , (2.15) 

does not give any simple criterion for complete positivity. Note, however, that ^'-representation 
is well suited for the composition of maps. If 



B = ^b',^E^^, (2.16) 

with 

b'^^ = {{B,E^^)) , (2.17) 
then the map C = Ao B gives rise to the following representation 

C = E^"/3^"/3' (2.18) 

where the matrix d = a' -h' . 
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Consider now a linear map A from C{Md, M^) and let us assume that A is diagonalizable, 
that is, it gives rise to the Jordan representation with 1-dimensional Jordan blocks. One has 

A = VDV-^, (2.19) 

where D is diagonal. It means that there exists an orthonormal basis fa in such that 

{{fa,Dfp)) =da5af, , (2.20) 

with da & C It shows that 

D=J2daPa, (2.21) 

Q!=0 

where 

Paa = fa{fa,a) , a £ . (2.22) 
Note, that a set Pa defines a family of orthogonal projectors 

PaPp = ^apPa , (2.23) 

together with 

^ Pa = id , (2.24) 

where id denotes an identity map in C{Mfi, M^). Hence, one obtains the following representation 
of ^ 

Aa = V DV-^a= ^ daVPaV'^a 

a=0 

= ^ d« VfaUa, V'^a) =Y,da V faiV'^* fa, O) . (2.25) 
a=0 Q!=0 

Let us define new basis 

ga := Vfa , ha := V-'*fa ■ (2.26) 
Note ga and ha define a pair of bi-orthogonal (or damping [19J) basis in 

{9a, hp) = {Vfa, V-^*fp) = {fa, f(s) = Sa^ • (2.27) 

Finally, one obtains the following spectral representation of the linear map A 

cP-l 

A=Y,daPa, (2.28) 

a=0 
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where 

Paa := ga{ha, a) , a e Ma . (2.29) 
Note, that a set Pa satsfies ^ ^ ^ 

PaPfi = SafiPa , (2-30) 

together with 

J2Pa = id. (2.31) 

a=0 

However, contrary to Pa operators Pa are not Hermitian, i.e. Pa 7^ Pa 

P*a:=K{ga,a) , a e . (2.32) 

It shows that Pq, are not projectors unless ga = ha- The corresponding spectral representation 
of the dual map reads as follows 

A#=^d;P#, (2.33) 

a=0 

where x stands for the complex conjugation of the complex number x. Hence, one obtains the 
following family of eigenvectors 

Aga = daga , A* ha = daha ■ (2.34) 

Consider for example a special case with V = U and J7 is a unitary operator in M^. One has 
U"^ = and hence = U. One obtains 

ga = ha = Ufa , (2.35) 

and hence P* = Pq. Note that 

A=Y,daPa, A*=Y, d^Pa , (2-36) 

q;=0 a=0 

which implies that the super-operator A is normal 

AA* = A* A . (2.37) 

3 How to generate commutative dynamics 

(k) 

Consider a family of Markovian semigroups A), defined by 

Af) = e*^'=, k = l,...,n, (3.1) 
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where are the corresponding generators. Suppose that are mutually commuting and define 

n 

A*o=Ep'^(*-*o)^So > (3-2) 
k=l 



where pfc(t) denotes time dependent probability distribution: pk{t) > andpi(t) + . . = 1. 

Let us observe that At^to defines a commutative non-Markovian evolution satisfying local in time 
Master Equation [18] 

At,to = ^t-tQAt,to , Ato,to = id . (3-3) 
To find the non-Markovian generator Ct let us assume the following spectral representation of 

LkP = Y.Xi'^gatT{h*,p) ■ (3.4) 

a 

One obtains 

CtP = ^ pa{t)gatTc{h*^p) , (3.5) 



with 



Hence, the solution At has the following form 



Pa{u)du] gatr{h*ap) . (3.7) 



Actually, one can easily generate a family of commuting generators Li, . . . Suppose one 

is given a Markovian generator L of a unital semigroup At = e^*. Denote by Ag the Laplace 
transform of At 



OO 1 
St ■ ^ 



As= / e-''Atdt = . (3.8) 

Jo s-L 

It is evident that for s > 0, ^4^ is completely positive. Moreover, 

:= sAs , (3.9) 

is unital. Indeed, one has 

/•OO 

$i°h = s / e-'^dtl = l . (3.10) 

JO 

Now, let us define 

fc+l jk 

One gets 
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(k) 

It is clear that for s > , is completely positive and unital. Therefore, for any integer k 
and s > one obtains the following Markovian generator 

L^'^) = $(*^) - id . (3.13) 

Hence, fixing s, one arrives at L/. := Ls^\ 

Let us observe that the construction of the commutative (s, A:)-family L^*^ may be used to 
construct a huge family of commuting time dependent generators. Note, that taking a discrete 
family of function 

/fc : M+ X M+ ^ M+ , 

one may define 

poo 

A[f]:=V/ fk{t,s)LfUs, (3.14) 
k 

where we used a compact notation f = (/i,/2,...). It is clear from the construction that 
[>Ct[f],>Cs[f]] = 0, and hence one easily find for the evolution 

At,t, [f] = exp ^* *° Cu[i] du^ . (3.15) 
defines a family of commuting time dependent Markovian generators. 

4 A class of commutative dynamics of stochastic classical sys- 
tems 

I Markovian classical dynamics 

Consider the dynamics of a stochastic d-level system described by a probability distribution 
P = (p(0)) • • • :P{d —I))- Its time evolution is defined by 

d-l 

Pt{m) = ^ Ttfy (m, n)po{n) , (4.1) 

n=0 

where (n, m) is a stochastic matrix satisfying the following time-dependent master equation 

Tt,to = LfTt^to ) ^to,to = , (4.2) 

that is 

d-l 

Tt,to{rn,n) = ^ Lt{m,k)Ttfy{k,n) , Tto,to{m,n) = S{n,m) . (4.3) 
k=o 

Let us assume that Lt defines a commuting family oi d x d matrices, i.e. 

d-l d-l 

^ Lt{m, k)Lu{k, n) = ^ L^im, k)Lt{k, n) , (4.4) 

A;=0 fc=0 
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for any t,u > to. A particular example of commutative dynamics is provided by circulant 
generators. Let us recall that a d x d matrix L{m,n) is circulant [20j if 

L(m, n) = a{m — n) , mod d , (4-5) 

that is L is defined in terms of a single vector a = (a(0), . . . , a{d — 1)). 

Proposition 1 Circulant matrices define a commutative subalgebra of M^. Hence, if L and L' 
are circulant then L" = LL' = L'L is circulant. Moreover, if 

Hhj) = a{i - j) , L'{i,j) = a{i - j) , L"{i,j) = a"{i - j) , mod d , 

then 

a" = a*a' , (4.6) 
where a* a' denotes a discrete convolution in Z^, i.e. 

d-l 

a"{n) = Y^a{n- k) a'{k) . (4.7) 

fc=0 

Therefore, multiplication of circulant matrices induces convolution of defining d-vectors. Inter- 
estingly, spectral properties of circulant matrices are governed by the following 

Proposition 2 The eigenvalues Im and eigenvectors ipm of a circulant matrix 

Llpm = Imlpm , (4.8) 

read as follows: 

d-l 

lm = Y.ak\^^ , (4.9) 

fc=0 

and 

{^^)^ = 1=X"^- , (4.10) 

where 

A = e^^^/^^ . (4.11) 

Let us observe that the Kolmogorov conditions for the stochastic circulant generator Lt give rise 
to the following condition upon the time dependent vector at{m): 

1. at{m) > 0, for m 7^ 

2. at(0)<0, 

3- E.m«t("^) = 0' 
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for t >to. Now, it is clear that the solution to (14.3 



Tt,to = exp (^^ Ludu^ , (4.12) 

defines a circulant stochastic matrix. Hence 

Ttfy {m, n) =: Pt^to {m - n) , (4.13) 
defines a time-dependent stochastic vector Pt^to i^) ■ Note that 

d-l d-l 



Pt{m) = '^Ttfy{m,n)po{n) = Pt,to{m - n)po{n) , (4.14) 

n=0 

Pt = Pt,to*Po- (4.15) 



n=0 n=0 

and hence 



One obtains from (j4.3p 

d-l 



^^b^ = Y^at{m-k)Pt,t,{k) , Pt,,to{m) = 6^0 , (4.16) 

which can be rewritten in terms of discrete convolution 

Pt,to =at* PtM ' PtaM = e , (4-17) 
where 'e' corresponds to the distribution concentrated at 0, i.e. e(m) = 5mO- 

Proposition 3 A convex set Vd of probabilistic d-vectors defines a semigroup with respect to 
the discrete convolution. The unit element e = (1,0, ... ,0) satisfies 

p*e = e*P = P , 

for all P £Vd. 

To solve ()4.17p one transform it via discrete Fourier transform to get 

dPf f,, (m) ~ , , , , 

'^^^ ' = at{m) Pt,t,{m) , ^toA.M = 1 , (4-18) 

where 

d-l 

x(n) = ^ A"'^ x(fe) , (4.19) 



A:=0 
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and the inverse transform reads 

d-l 

d 



x{k)='^^\-'^^x{n) . (4.20) 



n=0 

The sohition of (j4.18p reads as fohows 

(■t 



Pt,to{'m) = exp (^j^ auim)dt?j , (4.21) 



and hence one obtains for the stochastic vector Pt^to{fn) 

k=o ^•^*o 
It is clear that Pt^to satisfies the following composition law 

Pt,s * Ps,u = Pt,u , (4.23) 

or equivalently 

Pt,s ■ Ps,u = Pt,u , (4.24) 
for all t > s > u. In particular when a{n) does not depend on time then (I4.2ip simplifies to 

Pt,to (m) = exp {a{m) [t -to]) , (4.25) 

and hence 1-parameter semigroup Pt-to ■— Pt,to satisfies homogeneous composition law 

Pt*Ps = Pt+s , (4.26) 

or equivalently 

PfPs = Pt+s , (4.27) 

for all t > s > to- 
ll Non— Markovian classical dynamics 

Consider now the non-Markovian case governed by the following local in time master equation 

Pt,to = at-to * PtM > PtoM = e • (4-28) 
One easily obtain for the solution 

PtM M = IY1 (/ au{m)di^ . (4.29) 



11 



Note the crucial difference between (j4.22p and (I4.29p . The former defines inhomogeneous semi- 
group whereas the latter is homogeneous in time (depends upon the difference 'i — to') but does 
not define a semigroup, i.e. does not satisfy the composition law (I4.23p . 

Let us analyze conditions for which do guarantee that Pt^to defined in (|4.29|) is a probability 
vector, that is, 

d-l 

Pt,to{m)>0, ^Pt,to(m) = l, 

m=0 

for all t > tQ. It is clear from (j4.28p that a(T) has to satisfy 

au{m)du > , (4.30) 







for m > 0, and 

d-l 



y / au{m)du = , (4.31) 



m=0 ° 

which implies that 

/ au{^)du < , (4.32) 
Jo 

for all r > 0. These conditions generalize Kolmogorov conditions in the inhomogeneous Marko- 
vian case. We stress, that au(m) needs not be positive (for m > 0). One has a weaker condition 
(|4.30p . Note, that if at{in) > for m > 0, then au{m)du defines a monotonic function of 
time and hence the non-Markovian relaxation exp(Jp^ a„(m)(in) is monotonic in time as well. 
Finally, let us consider the corresponding nonlocal equation 

Pt^to = / Kt-u * Pu,to du , Pio,to = e , (4.33) 

J to 

with the memory kernel Kt-u- Note, that we already know solution represented by (j4.29p but 
still do not know the memory kernel K. Performing discrete Fourier transform one gets from 

Pt.to{'m)= Kt-u{m)Pu.to{m) du , Pt^^toim) = I ■ (4.34) 
Jto 

Define the time-dependent vector 

ft{m) = at{m)exp^J au{'m)dv^ , (4.35) 
then following [H] one obtains 

Ks{m) = -^iH- , (4.36) 
1 + fs[m) 

where Xg denote the Laplace transform of Xf. Clearly, the problem of performing the inverse 

Laplace transform Ks{m) — > Kt{m) is in general not feasible. Hence, the memory kernel 
remains unknown. Nevertheless, the solution is perfectly known. 
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Remark 2 Note that a stochastic map Po ^ Pt = ^t.toPo may be rewritten in a 'quantum 
fashion' as follows. Any probability distribution p = (p(0), . . . ,p{d — 1)) gives rise to a diagonal 
density matrix 

d-l 

P = ^p(n)e„„ , (4.37) 

n=0 

and the map po — 5- Pt reads as follows 

d-l 

Pt= ^ Tt,to{'m,n)emmPoenn ■ (4.38) 

m,,n=0 

III Dynamics of composite systems 

Consider now dynamics of A^-partite system living in = Tj^x . . . xZ^ . Let n = (ni, . . . , n^r), 
with n/j G Zrf and let 

Ptfy : ^ [0,1] , (4.39) 
be a probability vector living on satisfying the following Markovian master equation 

Pt,to = at * Pt,to , Pto,to=e, (4.40) 

where 'e' is defined by 

e(n) = (5no := 6^0 ■ ■ ■ 5nj^o ■ (4.41) 
Now, performing the discrete Fourier transform one gets 

^^^^1^ = 5t(m) Pi,i„(m) , Pto,to(m) = 1 , (4.42) 

where 

i(m) = ^A'"^rE(k) , (4.43) 

k 

and the inverse transform reads 

^(k) = J^E^^^'^^H • (4.44) 

m 

The solution of ()4.42p reads as follows 

Pt,io(m) = exp a„(m)dn^ , (4.45) 
and hence one obtains for the stochastic vector ^^^^^(m) 

P,,,„(m) = ^ E ^""'^ (yjyu{m)dv^ . (4.46) 
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It is clear that Pt^to satisfies the inhomogeneous composition law (I4.23p . If am is time independent 
then Pf^to defines 1-parameter semigroup Pr := Pr+to,to satisfying homogeneous composition law 

Km . ' 

Note, that in the case of non-Markovian dynamics one has 

Pt,to = at-to * PtM ' -fto.to = e > (4-47) 
giving rise to the following solution 

ft— to 



Pt,t, {^) = ^Y. ^""'^ ( / Sm(n)dn ) . (4.48) 



k -0 

The non-Markovian dynamics is time homogeneous but does not satisfy (j4.23p . 

5 A class of commutative quantum dynamics 

Consider now an abelian group x Z^. Equivalently, one may consider a cyclic toroidal lattice 
X Trf, where 

Trf = {A'" , m = 0,l,...,d-l} , (5.1) 

which is an abelian multiplicative group. Let us define the following representation of T,^ x 
in Md-. 

ZdxZdB (m, n) Umn e , (5.2) 

where Umn are unitary matrices defined as follows 

where {eo, . . . , e^-i} denotes an orthonormal basis in C^, and A stands for dth root of identity 
(see formula (14. IIP ). 

Proposition 4 Matrices Umn satisfy 

and the following orthogonality relations 

ix{u*^^uu) = ddmk^nl ■ (5-6) 

Hence, formula ()5.2p defines a projective representation of the abelian group x T^- It is 
therefore clear that 

Zrf X Zrf 9 (m,n) Umn G C{Md,Md) , (5.7) 

with 

Umna ■■= UmnaUmn , a e , (5.8) 
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defines the representation of Td x in the space of superoperators C{Md, Md). 
Now, for any 

a : Zd^^d C , (5.9) 

let us define a linear map A £ C{Md, M^) 

d-l 

A= a{m,n)Un,-m , (5.10) 

m,n=0 

that is, we define a representation of Aid in ^{Md, Md). 

Proposition 5 Ifa{m,n) G M, then A is self-adjoint, that is 

Ax* = {Ax)* , xeMd . (5.11) 

// a{m, n) > 0, then A is completely positive. If moreover ^ a{m, n) = 1, then A is trace 
preserving and unital. 

One proves the following 

Proposition 6 Let a,b,c £ Md be represented by A,B,C £ C{Md,Md), respectively, that is 

d-l d-l d-l 

A= a{m,n)Un-m , ^ = b{m,n)Un-m , = ^ c{m,n)Un-m ■ 

m,n=0 m,n=0 m,n=0 

Then A o B = C if and only if c = a*b. 

Hence, the set of maps constructed via (j5.1U|) defines a commutative subalgebra in C{Md, Md). 

Proposition 7 The spectral properties of the linear map 15. 1 0|) are characterized by 

Auki = OkiUki , (5.12) 
A*uli = akiuli , (5.13) 

and hence its spectral decomposition reads as follows 

d-l 

A = a{m,n)Pmn , (5.14) 

m,n=0 

where P„in a projector defined by 

Pmn X = — tr(li,^^x) , (5.15) 

for any x £ Md- 
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In particular, if a{m, n) is real, i.e. A is self-adjoint, then one has 



A-Uki = akiuki , (5.16) 
A*uli = mull. (5.17) 



Note, that the action of A upon the basis Cij is given by 



d-l 



'^(^ij= a(m,n)A"(* ^'^ ei_m,j-m ■ (5.18) 

m,n=0 

Hence, diagonal elements satisfy define an invariant subspace in 

d-l 

Aeii= ^ a{m,n)ei-rn,i-m ■ (5.19) 

m,n=0 

Let Pt^to : X — ^ [0, 1] satisfy the following inhomogeneous master equation 

Pt,to = at * Pt,to , Pto,to = e . (5.20) 
Now, following ()5.10p . let us define 

d-l 

Ato= Yl Pt,to{m,n)Un,-^m , (5.21) 

m,n=0 

and 

d-l 

^t= Y at{m,n)Un-m ■ (5.22) 

m,n=0 

Then, Proposition E] implies the following local master equation for the dynamical map At^to '■ 



Atfy = CtAt^to , Ato,to = id • (5.23) 
Note, that the time dependent Markovian generator may be rewrite as follows 

^tP = ^Y Mm,n)(^[Un-m,pUn,-m] + K,-mP, < - j) , (5.24) 

m,n 

where Yl'm,n^rnn ■= J2m,n^rnn " ^00- Hence, recalling that at{m,n) > for (m,n) / (0,0), 
the above formula provides the Lindblad form of Ct- The corresponding spectral representation 
of the generator reads as follows 

d-l 

Ct= Y at{m,n)Pmn- (5.25) 

m,n=0 
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Note, that due to aj(0, 0) = 0, one has Ct^d = 0. The corresponding solution of (|6.1ip is therefore 
given by 

^t,to= exp f / au{m,n)du] Pmn ■ (5.26) 

m,n=0 ^-^^O / 

If 

Pt^to satisfies non-Markovian classical master equation 

Pt,to = at-^to * PtM ' PtoM = e , (5.27) 
then the quantum dynamical map At^to satisfies non-Markovian equation 

■At,to = ^t-toAt,to ) ^to,to — id , (5.28) 
with the solution given by the following formula 

Ato = / au{m,n)du] Pmn- (5.29) 

m,n=0 ^-^0 ^ 

This spectral representation of Ar := AtQ+T,to enables one to construct the corresponding mem- 
ory kernel /C,-. Using the following representation [15] 



Ar = id + / Fsds , (5.30) 
Jo 

where 

Fs = CsAs , (5.31) 
one finds the spectral representation for the super-operator function Fs'. 

Ft = ^Mm,n)Pmn , (5.32) 
m,n 

with 

ft 



ft{m,n) = at{m.,n) exp au{m,n)duj . (5.33) 
Therefore, one may write the corresponding non-local equation 

At= f ICt-uAudu , (5.34) 
Jo 

with the memory kernel is defined in terms of its Laplace transform as follows 



p \ ^ sfs{m,n) 

iCs — / ^ ~ Pmn ) (5.35) 

m,n l + /s(m,n) 

where fs{Tn,n) denotes the Laplace transform of ft{m,n). Note, that in general one is not able 
to invert the Laplace transform ICs and hence the above formula in general does not have any 
practical meaning. 
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6 Dynamics of composite quantum systems 

Consider now a quantum dynamics of A^-partite d-level quantum systems defined by 

A,to= Yl Pt,to{m,^)Un-m , (6.1) 

where 

Uk,\x = Uk^ixul i , (6.2) 

for xeMf^, and 

Uk,\ = Uki,h <8) . . . (8) Ukt,,i,^ . (6.3) 

Proposition 8 Matrices Um,n satisfy 

^m,n^^r,s ~ ^m+r,n+s j (6.4) 

^m,„ = A™" (6.5) 
and the following orthogonality relations 

tr(^^m,n^k,l) = 6m,kSn,\ ■ (6.6) 

The spectral representation of At^t^ has the fohowing form 

At,to= J2 ^t,to(m,n)P„,„ , (6.7) 



where Pm,n is a projector defined by 

nr* — 



Pm,nX — jAjT '^m,n tr('Uni nX) , (6.8) 



for any x G . Assuming that 

Ptfy, : X ^ [0, 1] , (6.9) 
satisfies classical Markovian inhomogeneous master equation 

PtM = * PtM > PtoM = e , (6.10) 

one obtains 

At^tQ = ^tAt,to , Ata,tQ = id , (6.11) 
where the time dependent Markovian generator is defined by 

^t= Yl 2t(m,n)P™,„ . (6.12) 

Hence, the corresponding solution reads as follows 

At,to= Y ^^P f / a«(m,n)du'j Pm,„ . (6.13) 
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7 Commutative dynamics of 2-level system 

Consider the time dependent generator for a 2-level system defined by 



= -^e(i)k3,P] +7(i)(/^^l + (1 - /^)'C2)/3+ ^ J2 Cal3{t)[[TTa,p7ri3] + [TTaP,T^^]) , (7.1) 

a,/3=0 

where the time independent Markovian generators £i and £2 are defined as follows 

C2P = a'pa^ -^{a^a',p} . 



One easily shows that 

[A,£s]=0, (7.2) 

and hence Ct does generate a commutative quantum dynamics. In (17. ID the 'mixing' parameter 
/i G [0, 1], and projectors tTq are defined by 

ttq = a^a^ , TTi = a^a~ . (7.3) 

Note, that if 7(i) > and the time dependent matrix ||co^(t)|| is semi-positive definite, than Ct 
defines time dependent Markovian generator. If 

t 

-i{u)du > , (7.4) 





and the matrix 

ft 



Cai3{u)du 1 1 > , (7.5) 



for all t > 0, then Ct generates non-Markovian dynamics. 

One easily solves the corresponding spectral problem for Ct 

Ctoj = , 

Ct a+ = T{t) a+ , 

CtCF^ = r(t)fT" , 

A 0-3 = -i{t) 0-3 , 

where the invariant state u reads as follows 

uj = piTi + il- p)tto , (7.6) 
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and the complex eigenvalue r{t) is defined by 



r(0 = -i l{t) + coo{t) + cu{t)-2cio{t) + 2is{t) . (7.7) 



2 

Similarly, one solves for the dual generator 

4h = 0, 

Ji:fa+ = ma+ , 

jrfa- = r{t)a- , 

jrfa = -^{t)a, 

where ^ 

(7 = (1 - /x)7ri - /XTTo = 2 (^^3 - Il2tr(a;(73)^ . (7.8) 

Hence, introducing a bi-orthogonal basis 



such that 
one has 

with 



90 


= OJ 


ho 


= h 


91 


= 0-"*", hi 


= a' 


92 


= a 


h2 


= a 


93 


= 0-3 


hs 


= a 


{9a, 


h0) = 


^^i9ahp) 


= Sat 



(7.9) 

3 

ji:tP = Y.Xa{t)9atT{hlp) , (7.10) 



a=0 



Ao(t) = 0, Ai(i) = A2(t) = r(t) , A3(t) = -7(i)- (7.11) 
Hence, the solution to the Markovian master equation 

At^to = A^t,to , Ato,to = id , (7.12) 
reads ^ 

AtoP = y]exp(/ Xa{u)du] gatT{hlp) . (7.13) 

a=0 ^-^^0 J 

Consider now 

V : M2 — ^ M2 , (7.14) 

defined by 

V a = eoo (^fJ- tr(eiia) + tr(eooo)) + en (^{1 - fi) tr (eno) - tr(eooa)) 

+ eiotr(eoia) + eoitr(eioa) . (7.15) 
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One easily finds for the inverse 

y-^a = eoo( - Aitr(eiia) + (1 - /i)tr(eooa)) + eii(^tr(eiia) +tr(eooa) 
+ eiotr(eoia) + eoitr(eioa) , 

and hence 

V~'^*a = eoo(tr(eiia) + (1 - /i)tr(eooa)) + eii(^tr(eiia) - /itr(eooa) 
+ eiotr(eoia) + eoitr(eioa) • 

One finds 

^600 = 0-3, Veii=uj, V = , 

and 

y-i#eoo = a , V-^*eii = I2 , V-^*a^ = . 

Hence, defining 

/o = en , /i = 0-+ , /2 = fj" , /s = eoo , 

one has 

ga = V fa , K = V~^*fa , 

which shows that V diagonahzes Ct and ^t,to, that is, 

3 

Ct = Y,K{t)vPay-^ , 

and 

3 



^— n \-Jto J 



where 

P.P = fMflp) 

8 Conclusions 



7.16) 

7.17) 

7.18) 

7.19) 
7.20) 
7.21) 

7.22) 

7.23) 
7.24) 



In this paper we analyzed a class of commutative dynamics of quantum open systems. It is 
shown that such evolution may be effectively described via spectral analysis of the correspond- 
ing time dependent generators. The characteristic feature of the corresponding time-dependent 
dynamical map is that all its eigenvectors do not depend on time (only its eigenvalues do). 
Actually, majority of examples studied in the literature (see e.g. [I]) belong to this class. If 
eigenvectors vary in time then the solution is formally defined by the time ordered exponen- 
tial but the problem of finding an explicit solution is rather untractable. We stress that both 
Markovian and non-Markovian dynamics were studied. Our analysis shows that the local ap- 
proach to non-Markovian dynamics proposed in |18j is much more suitable in practice than the 
corresponding non-local approach based on the memory kernel. 
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